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Abstract 

We investigate the relevance of the metric and of the geometry in five-dimensional 
models of hadrons. Generically, the metric does not affect strongly the results and even 
flat space agrees reasonably well with the data. Nevertheless, we observe a preference for 
a decreasing warp factor, for example AdS space. The Sakai-Sugimoto model reduces to 
one of these models and the level of agreement is similar to the one of flat space. We 
also consider the discrete version of the five-dimensional models, obtained by dimensional 
deconstruction. We find that essentially all the relevant features of "holographic"' models 
of QCD can be reproduced with a simple 3-site model describing only the states below 
the cut-off of the theory. 
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1 The Road to Five Dimensions 



The search for a dual description of QCD is a task of formidable difficulty and importance: a 
success would greatly improve the still limited understanding of qiiantum field theory faraway 
from the perturbative regime. Such a dual model, provided it exists, should be weakly coupled 
in the IR where QCD is strongly coupled and the relevant degrees of freedom are hadrons. In 
the IR, therefore, the dual model should reduce to a calculable effective (field or most likely 
string) theory of hadrons and one might try, guided by the observations, to guess its broad 
features following a bottom-up approach. 

Apart from the far IR, where the only relevant degrees of freedom are the Goldstone bosons 
whose dynamics is entirely determined by the QCD global symmetries, the very existence of such 
a model can be doubtful. Indeed one observes no sharp separation between the mass and size 
scales of the hadrons, which seem to be all fixed by the dynamically generated confinement scale 
of QCD, Aqcd ~ 1 GeV. The problem is that no obvious weak couphng can be identified from 
the observational point of view, while a weak coupling is the starting point for the construction 
of any effective theory. 

This difficulty immediately shows up in concrete attempts of describing hadrons, and in 
particular the lightest spin-one meson, the p(770). Both in the popular and phenomenologically 
successful Hidden Local Symmetry (HLS) approach [1] (see [2] for a recent review) and in the 
less successful (but theoretically compelling) Georgi's "vector" model [3], agreement with data 
requires a near to maximal coupling {gp ~ 6) for the p meson. Even though a perturbative 
expansion can be set up and loop corrections computed [2] by treating gp formally as a small 
expansion parameter, the fact that gp is numerically so close to the perturbative bound of 
Air makes us doubt that the entire approach makes sense. Is the fact that gp < Att simply 
an accident or is there a parametric reason why this happens? More precisely: is there a 
deformation of the QCD theory in which gp ^ and treating gp ets a weak coupling is justified? 

The answer to both questions is affirmative and has been found in the 't Hooft large-A^c 
limit [4]. At large- A'c, assuming that confinement persists, a weak coupling emerges (see [5] and 
references therein) and QCD becomes a weakly coupled theory of mesons. The meson masses 
have a finite large- A^c limit while the couplings scale as 1/a/A^; this suggests that a weakly 
coupled effective description of large-Ac mesons must exist. It is hard to believe, however, 
that such a description is provided by the models mentioned in the previous paragraph. They 
describe indeed a single resonance while it is known that large- Ac mesons necessarily arise in 
infinite towers [5]. For those models to make sense at large- A^c one should imagine, as suggested 
in [3], the higher resonances to be much heavier than than the first one, or almost decoupled. 
This seems unlikely to happen. 



1 



In this context, the recently proposed 5d models of mesons [6-9] (see also [10] for previous 
attempts) represent a clear progress. In these models, usually denoted as "holographic QCD" or 
"AdS / QCD" , the global chiral symmetry of QCD is promoted to a local gauge symmetry in a 5d 
bulk and its spontaneous breaking is implemented by either a 5d Higgs mechanism [7-9] or by 
boundary conditions [6,11]. Both cases, which actually give very similar results, automatically 
deliver infinite towers of 4d vector meson fields of fixed isospin quantum numbers. Not only this 
structure puts us closer to the large-iVc expectations, it also makes the model more predictive: 
the "minimal" scenario with breaking by boundary conditions only requires 3 parameters at 
tree-level and describes at the same time the p(770), cij(782), ai(1260) and /i(1285) mesons 
and their towers, plus of course the pions. [] With the same set of parameters, the minimal 
model also describes baryons which, as expected at large-Ai'c, are described as calculable 5d 
skyrmions [12,13]. 

Another important feature of the 5d construction is that it provides a simple bookkeping 
of the Nc factors: the Nc scalings of the various observables are recovered, both in the mesonic 
and baryonic sector, if taking the 5d expansion parameter {i.e. the gauge coupling g^) to scale 
as (^5 oc 1/ v^iVc. Notice however that the quantitative agreement of the model with observation 
(which can be quantified as 10% in the meson sector and 30% in the baryon one) cannot be 
ascribed to the combined use of naive laige-Nc considerations and of chiral symmetry, the 5d 
model gives much more precise informations. It implies for instance certain sum rules [6,11] 
that we will also discuss in the following and, remarkably, makes Vector Meson Dominance 
(VMD) arise automatically. 

The absence of higher spin states, and of the entire Regge phenomenology, is the main 
limitation of the 5d approach (see however [15] for a discussion of this and related topics), even 
though this does not result in a concrete phenomenological failure of the model. For the real- 
world case of Nc = 3, indeed, the regime of validity of the 5d model is limited and the observed 
higher spin states, which are quite heavy, could be thought to be above the 5d cutoff. We 
should not worry, for the same reason, about the masses of the heavy resonances in each tower, 
which will quickly become too broad to be identified as particles and merge in the continuum; 
ordinary NDA considerations suggest that it should be difficult to go above the p(1450). At 
the theoretical level, the problem shows up for larger values of A^^^ where the coupling decreases 
and the cutoff grows, but still higher spin states do not appear. For the model to make sense 
as a description of large- A'c we should assume higher spins to be parametrically heavier, which 
seems difficult, or mostly decoupled from the vector mesons. This seem unlikely but, according 
to what we know, is not in contrast with any large- argument. 

The phenomenological road to five dimensions leaves the metric of the 5d space completely 

^In non-minimal models, scalar resonances are also present [9]. 
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unspecified. While originally inspired by the AdS/CFT correspondence [16], the discussion 
above makes it clear that the AdS/QCD approach can be motivated by bottom-up and phe- 
nomenological considerations which do not require AdS geometry. Indeed, especially at zero 
temperature, it appears difficult to use AdS/CFT for the strong interactions: QCD is conformal 
at weak coupling while the duality becomes useful in the opposite regime, and it is not con- 
formal at large distances, i.e. in the regime relevant to determine the properties of hadronsj^ 
All this motivates us to consider five dimensional models of hadrons with general warp factor. 
As long as the 4d Poincare invariance is preserved, any metric would lead to a "reasonable" 
5d model, enjoying the general features mentioned above. For example, the popular Sakai- 
Sugimoto model [18] is equivalent, for what concerns the physics in the meson sector, to an 
AdS/QCD model of this kind with a geometry different from AdS. 

Since certain relations among observables will be metric-independent, all geometries will 
share certain common predictions, but the agreement with data of different models will ulti- 
mately depend on the shape of the metric. In the rest of the paper we investigate quantitatively 
this dependence. The main result is that, even though "very wrong" choices can be done, a 
large class of metrics, including flat space, produces a good agreement with the data. At first 
sight, considered the large expected corrections to our leading order results, no metric among 
fiat space, AdS or Sakai-Sugimoto appears to be strongly favored. Taking the error more se- 
riously, AdS seem to be preferred, but its phenomenological success is shared by many other 
metrics, for example a warp factor ^JL/ z. Generically, a mildly decreasing warp factor provides 
a good fit of the data. Another important feature is that the 5d models have an "intrinsic" 
minimal error which depends on the successful, but not perfect, model independent predictions. 
The error is almost minimized by the AdS or \/ Ljz geometries. 

Having quantified the dependence on the metric of the AdS/QCD results, the natural ques- 
tion is whether the role of the fifth dimension, i.e. of the 5d geometry, is itself crucial. In 
fact, since the 5d theories required to reproduce QCD are effective field theories with a limited 
regime of validity, we might expect that a model with only the few modes below the cut-off 
(which turns out to be of the order 3-4) will be sufficient to reproduce similar results. Such 
a model can be obtained by deconstructing the ffith dimension [6, 19]. In the deconstructed 
version the five dimensional gauge theory is replaced by a 4d gauge theory with K sites and 
nearest neighbor interactions. We will consider models with different number of sites and show 
that already the case of if = 3 provides an good fit of the data (better than fiat space) while 
if = 4 is as precise as AdS. Moreover, since in the limit if — > oo the discrete model is classically 
equivalent to the 5d geometry, we use the deconstructed theory as a tool to numerically solve 
arbitrarily complicated metrics. 

^See however [17] for the opposite claim. 
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The paper is organized as follows. In next section we review the five- dimensional theory 
relevant to describe the spin 1 mesons of QCD and its deconstructed version. In section 3 
we present the results of the fit for different low energy QCD observables for continuous and 
discrete models. The validity of these effective theories is discussed in section 4. We conclude 
in section 5. The appendix contains a detailed derivation of the analog of the Chern-Simons 
term in the discrete models which is necessary to reproduce the anomalies. 



2 Models of Hadrons in 5d and 4d 
2.1 Continuous Models 

We consider the same model as in [13, 14], a pure five-dimensional Yang-Mills theory with 
U{2)lxU {2)r gauge group and warped extra dimension z G [z^v, z\r\. We denote as L = Zi^ — z-^v 
the "conformal length" of the extra dimension and by Ijm and Rm the U{2)l and U{2)ff gauge 
connections. The metric on the 5d space is 

gMN = a{zy r]MN , (2.1) 

where M and run over the five space-time indices, tjmn is the 5d Minkowskian metric with 
mostly minus signature and a{z), the warp factor, is a positive and regular function. Indexes 
are raised and lowered with the fiat metric. Any constant rescaling of the warp factor can 
be reabsorbed into the definition of the 5d gauge coupling g^, so that we can normalize it to 
a{z[^) = 1. This chiral gauge symmetry is broken at the IR boundary by the conditions 

(L/x — R/x) \z=zi, = , (L^5 + R^s) = , (2.2) 

where the 5d field strength is Lmn = "^a/La? — djyLM — ^[Lm, Ljv] and analogously for Hmn- 
At the UV boundary the fields can be identified with the external sources for the QCD global 
currents 

^^l\z=zuv ~ -^M ' R^t |2;=2uv = I'm- (2-3) 

We will consider general functional forms of a{z), including for instance fiat space, a{z) = 1. 
The geometry, therefore, does not distinguish between the UV and IR boundaries, and what 
makes them different are the boundary conditions ( 2.2|2.3 ). 



Thanks to the IR boundary conditions (2.2) we can reformulate our model, which describes 



two f/(2) gauge fields L and R on the space z G [z^^, ^ir], in terms of a single U{2) field A living 
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on a "doubled" 5d space with extra coordinate u G [—L,L]. This is done by "gluing" L and R 
together in the following way: 



-R-/^ (^ir 



u e [-L,o] 

u G [0, L] 

L^^Zir + u), u G [-L,Q] 
-R5(2;ir - m), m G [0,L] 



(2.4) 



Notice that the IR boundary conditions (2.2), which L and R have to respect, ensure that the 
field A is continuous at m = 0. 

The doubled space has two boundaries u = —L and u = L which we call respectively the 



left and right boundaries. The boundary conditions for A, obtained by combining eq.s (2.3) 

A^(-L) = 1^, A^(L) 



and (2.4), are 



(2.5) 



The metric of the doubled space is of the form (2.1) with a warp factor 



a[u] 



a{zi^. + n), n G [— L, 0] 
a^Zij. — u), M G [0, L] 



(2.6) 



which is symmetric under the reflection u —u. Being this transformation an isometry, we 
could impose it as a symmetry of our model which however would not correspond to any of the 
QCD symmetries. |^ What corresponds to the QCD parity is, in this language, the simultaneous 
inversion of all spacial coordinates, i.e. the combined action of u — > —u and x — > —x where 
X denotes 3-space coordinates. In the language of L and R flelds this symmetry corresponds, 
via eq. (2.4), to the L ^ R interchange combined with x 



Finally, the action reads 
1 



YM 



d^x a{u) 



a 

y 



—X. 



cs 



167r2 



d^x 



MNOPQ 2 P P 



(2.7) 
(2.8) 



where the fleld strength is Fmn = QmAn — On Am — i[AM,AN], we have parametrized the 
U{2) gauge fleld A as A^ = + Am^/2 in terms of Pauli matrices and (.) denotes 

the trace. In the previous action we introduced different couplings for U{1) and SU{2) sectors 
as this is allowed by the symmetries of the model. However it is consistent with large Nc to 



^This would rather correspond to the U ^ W "accidental" symmetry of the ordinary pion Lagrangian, 
which is broken by the Wess-Zumino-Witten term. Similarly (see below) u — > —u is broken by the CS term. 
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set them equal, so we will work with a = 1 in the rest of the paper for convenience, since it 
would not affect the results much. The presence of the Chern-Simons term Scs is required for 
matching the anomalies of large-A^^c QCD and its coefficient is fixed by the number of colors 
Nc- For a given warp factor a{z), therefore, the model has 2 free parameters to be determined 
by experiments. |^ Notice that the simplified form of the CS that we use is only valid in the two 
flavors case, i.e. for U{2) 5d gauge fields. It is very simple, splitting the u domain of integration 



into u G [—L, 0] and u G [0, L] and using eq. (2.4), to show that the action (2.8) is equal to the 



one considered in [13] once expressed in terms of L and R fields. The two theories are therefore 
equivalent up to the choice of the warp factor which was fixed to a{z) = L/z in [13] while we 
will consider more general possibilities. 

One of the advantages of rewriting the model on the doubled space is that the connection 
with the Sakai-Sugimoto model becomes manifest. This model is motivated by strings but 
effectively reduces, in the limit of large string scale, to a gauge field theory on warped 5d space 
of the kind considered in this paper. All the calculations of mesons properties [18] are performed 
in the field theory regime, and are insensitive to the string completion. The action is 

Sym = [ d'xdz {F,,Fn + 2K{z) {F,,Fn] , (2.9) 



with K = 1 + z"^ plus the same Chern-Simons term as in (2.8). By the change of coordinates 



+ (2.10) 

(and rescaling u suitably) the action above can be rewritten in our frame with warp factor 

a{u) = K{z{u)y^ . (2.11) 

Therefore the Sakai-Sugimoto model is at the practical level equivalent to an AdS/QCD model 
of the kind we consider here. 

It would be straightforward, for a generic warp factor, to study the phenomenological con- 
sequences of the setup described in this section. It is only for few particular choices of the 
warp factor, however, that the model could be solved analytically. For a generic warp factor 
we should in any case rely on numerical methods to solve the differential equations which de- 
termine the masses and wave functions of the KK resonances. This motivates us to study the 
continuous 5d models as a limiting case of the discrete model to which we now turn. The prob- 
lem of finding wave functions and masses will be reduced to the one of diagonalizing "large" 



AdS/QCD is often regarded as a 1-parameter model since the coupling is "matched" at the UV to the 
perturbative QCD result [7, 8] . In the spirit of the introduction we instead consider as a purely IR parameter 
and regard as a mere coincidence the fact that its matched value would be close, in the case of AdS, to the 
phenomenologically preferred one. 
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matrices, a numerically easier task. The discrete model, however, is not just a technical tool. 
It can be considered per se, for a limited number of sites, as a model of hadrons. It is only for 
an infinite number of sites that it (classically) coincides with the 5d model. 



2.2 Discrete Models 



What we denote as discrete models can be considered to be generalizations of Georgi's "vector" 
model [3] to the case of many sites, i.e. of many non-linearly realized gauge groups. The 
models we are considering are therefore the ones discussed in ref. [6], but with U{2) instead 
of SU{2) groups and the addition of a gauged Wess-Zumino-Witten (WZW) term at each 
link. The latter terms are necessary to reproduce the QCD anomaly and provide, expectedly, 
a discretized version of the 5d CS term [20]. 

A diagrammatic representation of the discrete model is shown in figure [1} each of the K + 1 
links represents a U{2) cr-model matrix (A; = 1, . . . , + 1) transforming as 



t 

•k ' 



(2.12) 



under a {U{2)l x [7(2)^)^ global group of elements L^. and R^. The associated non-dynamical 
sources are denoted as 1'^ and r'^ in figure [T] For the time being, different a-models do not 
interact with each other and at two derivatives order the Lagrangian allowed by the symmetries 



(2.12) is 



K+l 2 



k=l 



(2.13) 



where we have separated each U{2) matrix S'' = S'^S'' in its SU{2) (S'') and f/(l) (E* 



components. As in eq. (2.8) we choose the parameters of the U{1) and SU{2) sectors to be 



equal. Notice that we could have enlarged the symmetry group of the theory if we had chosen all 



the fk to be equal in eq. (2.13). With this choice we would have gained the discrete symmetry 



of links permutation and greatly reduced the number of parameters. 

Let us now make different links interact by performing a weak gauging of a subgroup of 
the {U(2)l X U{2)nf^^ global group; this partially breaks the global symmetry which forced 
us to write the Lagrangian (2.13), but in a controlled way, as explained in [3]. The groups we 



want to gauge are the vector combination of the Right group of each link (z.e, in eq. (2.12)) 
with the Left group of the following one (z.e, L^+i). The gauging procedure is depicted in the 
lower part of figure [l| each gauged group is denoted by a circle and the associated gauge fields 
by Al. The gauging consists in introducing the covariant derivatives 



(2.14) 
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si 



1 = 11 ri P 



„K lif+i 



1 = AO 





r = 



Figure 1: A graphical representation of the discrete modeL 



and adding a kinetic term to the A^, the Lagrangian reads 

^ 1 

-T- 



K+1 r2 



k=l 



fc|2 



k=l 



(2.15) 



The QCD Left and Right chiral transformations are identified, respectively, with Li and 



Rx+i in eq. (2.12), which are left ungauged. The fields A° and A^+^ appearing in the covariant 



derivative of, respectively, 5]° and are defined as the Left (A° = 1) and Right (A^+^ = r) 

QCD sources and are therefore non-dynamical field. They couple to the currents. 



lu=0 



J^^^ = z%((D,S^+V(S^+^))" , (2.16) 



which we identify with the SU{2) QCD currents and analogous expressions for the U{1) cur- 
rents. 

In this model QCD anomalies can be reproduced by introducing WZW terms at each site 
[20]. We postpone the full derivation of this term to the appendix. The final result is simply. 



'967r2 



r 

/ J2 ^'''^ (^^''"^ + d^'') - ((dA^-^)2 + {dA^y + dA''-^ dl*^) 
k=i I 



3 



+(A 



A" 



(2.17) 



where II'^ is defined by S'^ = exp i U'^ and we used a form notation in which A = —iA^dx^ and 
F = dA + A^. Notice the addition of 5*^ to the action makes an extra contribution to the U{1) 
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currents (2.16) arise. The latter can however be neglected since it will not contribute to the 
observables we will be interested in. 

This completes the construction of the discrete model. We still need to impose invariance 
under the QCD parity operation, which we identify with the reflection with respect to the middle 
point of the diagram in figure[l| i.e. as Yi^ ^ (S^+^"^)^ and A'^ ^ p^K+i-k combined of course 
with ordinary x — > — x. The reason for this identification is that parity must interchange the 
Right source r = A^+^ with the Left one 1 = A°. In order for the action to be invariant we 
have to take fk = fK+2-k and gk = gx+i-k- 

It is important to remark that the logic we have followed in the formulation of the model 
led us to write a Lagrangian of special "nearest-neighbor" form 

C=^Ck (S^ A^-\ A^) . (2.18) 

k=l 

Non-nearest neighbor interactions will be radiatively generated by the gauge fields dynam- 
ics which breaks the symmetry (2.12) and makes different links interact. Their coefficients, 
however, will be suppressed by powers of the gauge couplings Qk, which are the parameters con- 
trolling the breaking of the symmetry (2.12). In an expansion in the gauge couplings, therefore, 
non-nearest neighbor operators are subleading and we are allowed not to consider them at the 
leading order. The situation is exactly the same as in the model of ref. [3] which corresponds 
to the K = 1 case. It is clear from the discussion above that the discrete model is based on a 
weak coupling (small gk) expansion, and that it is not meaningful to consider gauge couplings 
which exceed Att. We will however use formally infinite values of the couplings when we will 
take the "continuous limit" of the discrete model: we have already mentioned that this limiting 
procedure can only be thought as a numerically convenient method to obtain predictions in the 
5d theory with general metric. 

Following [6,19], the discrete model becomes equivalent to the continuous one described in 
the previous section when the the number of sites K is taken to infinity while keeping gu^ /K 
and fk' /K fixed. The relation between 4d and 5d models is given by 

2 _ j^K + l 2 _ a[{uk-i + uu)/2]2{K + l) 

~ aiuk) 2L ' ^'^ - g,^ L ' ^''''^ 

It is already clear that we will identify the parameters L and g^ and the function a{u) with 
the ones appearing in the continuos 5d model. Moreover, the Uk are points on the discretized 
extra-dimension u, whose values are 

Uk = Uo + -^^k, (2.20) 
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with k = 0, . . . , K + 1 and uq = —L. 

In the continuous hmit, the fields S'^ are identified with the hnk variables of lattice gauge 
theories, i.e. with the Wilson lines on a straight path from the Uk to the Uk-i points, while 
is the 5d gauge field at the point Uk- It follows that, in this limit 



11 - z 



2L 



— yA5[K_i + Ufc)/2] , 



A/, [uk] 



2L 



(2.21) 



(Aj-Aj-i) ^ (95A,)[K„i+«,)/2] . 

It is simple to check, using the equations above and changing the sum into an integral, that the 



Lagrangian (2.15) reduces to the 5d Sym action in eq. (2.7). Similarly Sa reconstructs, the 5d 



CS action of eq. (2.8). The first line of eq. 2.17 reduces, using eq. (2.21), to the second term 
of eq. (2.8). One has to use that A^^^A^ = [A'^^^ — A^)A^ and neglect terms which go to zero 
faster than 1/K. For the second line, since (S^)^i5^S^ ~ ~i2L/{K + 1)-F^5 one obtains the 



first term of eq. (2.8). Moreover, the third line vanishes in the continuous limit as it scales like 
l/K\ 



2.3 QCD Observables 

We now present the expressions for the QCD observables to be compared with the data. We 
will derive the formulae in the discrete model; the analogous expressions in the continuous limit 



can be quickly obtained by taking K ^ oo and using eqs. (2.19). 



2.3.1 Chiral Lagrangian 

The effective Lagrangian for the pions up to fourth order in the derivatives is customarily 
parametrized as follows [21] 

J TT 



L 2 



4 



Li ((D^S)^(D^S)) + L2 ((D^S)^(Z},E)) ((D'^S)^(D'^S) 
+ Ls ((D^S)t(D^S)(D,S)t(D'^S)) 



(2.22) 
(2.23) 



■^What had to be shown is that the only field configurations relevant in the continuous limit are in the form 
of eq. (2.21 ) where, for instance, S'' only deviates from the identity by an 0{\/K) term. This could be checked 



explicitely to hold for the eigenmodes that we will derive in the following section. 
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where the covariant derivative is -D^S = — i /^S + zSr^. The one above is the purely 
"5't/(2)" part of the Goldstone Bosons Lagrangian, and it will not receive contributions from 



the Sa term in eq. (2.17). We therefore ignore the latter in the discussion which follows. 

To extract the coefficients of the pion lagrangian from our action it is convenient to follow 
the "holographic" approach of ref. [22], appropriately adapted to the discrete case. To do this 
we first choose a gauge where all the link variables S'^ but the first one, denoted as S, are 
set to the identity and then notice that the dependence on S of the gauge-fixed action can be 



reabsorbed in the redefinition of A^^. Namely, the action can be written as in eq. (2.15) with 
all the S'^'s set to the identity but with 

Aj = St(/, + z9,)S, <+^ = r,. (2.24) 

In terms of the original fields (before fixing the gauge), S is 

S = S^S^.-.S^+i, (2.25) 

and represents the discretized version of the Wilson line from one boundary to the other. 

We will now write the effective theory for the S field, which describes the physical pion, 
by integrating out the gauge fields A'^ at tree-level, i.e. by solving the classical equations of 
motion and plugging back into the action. It is possible to check that, like in the continuous 
case of [22], it is enough to consider the linearized equation at zero momentum in order to 
obtain the effective Lagrangian up to 0(j)^). The soluton is 

4 = ^ (1 + ^o) ^Hl, + ^d,)E + 1 (1 - 69 r, , (2.26) 



where 



"5 = 1-215:^1 itjf]- (227) 




The 0{p ) and 0{p ) terms in the chiral Lagrangian arise, respectively, from the first and 



second term of the Lagrangian (2.15). We immediatly read the pion decay constant 

K+l 



= E 72 ' (2-28) 

•I TT J k 



k=l 

which in the continuos limit this reduces to the standard, 



2 f du 
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For the order terms one finds, 



k=l 



K 

I . 

Lg = —Lio 



iEi(l-(^of) ' (2-30) 

^ k=l 3k ^ ^ 



which again reduce to the standard 5d formulae in the continuous limit. 
2.3.2 Vector Mesons 

The other observables are more easily obtained performing the analog of the Kaluza-Klein 
reduction. To do this it is convenient to proceed as in ref . [6] . By choosing the gauge 



exp 



■k 
Jk 



k J 



(2.31) 



the scalars do not mix with the vector fields. We can then diagonalize the mass matrix for the 
massive gauge bosons. To do this we introduce 



K 



Ki^)=9kY.^'nB;{x), (2.32) 



n=l 



where B'^{x) is a vector field with mass m^. 



Masses and wave functions are determined by the following eigenvalue problem: 

9k {fk 9k-i b'-' - Ul + fl+i) 9k bl, + fl^, g^+i + 4m^ 6^ = , (2.33) 

for k G [1, K] and where 6° = b^^^ = 0, since the I and r sources are non-dynamical and have 
to be put to zero. 

Due to the symmetric nature of our action the modes split into even (6^ = b^^^~^) and 
odd (6^ = —b^^^^^) whose masses alternate in the spectrum. The actual solutions (with unit 
norm) can be easily found numerically for given fk and gt (in particular if we fix them with 



the relation (2.19) for a given metric). 



The trilinear vertices we are interested in are 









c - 






^B^n-y 


= -gB-.^-e^"'"' d^BlOpB^Ti^ 




c - 


= -9b^., ■ d,B2 d,B: . 


(2.34) 
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The first two couplings can be obtained by simply plugging the wave-functions (2.33) into 
the action. To obtain the last two coupling we need to gauge the electromagnetism. In the 
language of the previous section the electromagnetic current corresponds to the combination 
Jem ~ V3 Jy + Jy^ ■ To make the photon dynamical (see [13] for the discussion in the 
continuous case) we just need to add a kinetic term for the source associated to the current, i.e. 
for the appropriate components of Aq and Ak+i- We could now repeat the full Kaluza-Klein 
reduction. However it is much easier to recognize that to leading order in the electric charge 
the solution is simply 

K 



n=l 
K 



Alix) = \eBl{x) + Y, 9khlB;{x), (2.35) 



3 



which follows from the fact that the heavy modes are unchanged to leading order in the per- 
turbation. Using these wave-functions we obtain 

J 2 ^^+1 ^ _ , 



K+l 



9b^b^^. = ^ J2 {dk-i bt ' + gk bt^ [gk^i bl ^ + gk 6') , 



32-^ tt ft 

eN, 

24^/, 
eN, 



gB'^TT'j — — gs'^-KiT ^ 



9B"nf - g^2f 9B"777r ■ (2.36) 



Notice that the second coupling is non zero only if 6^ and 6^ have the same parity while the 
others only if is even under the refiection k K + 1 — k {i.e. for vector resonances). 

The last observables that we consider are decay constants (for axial and vectors mesons) 
and the pion radius. As for the meson decay constants are defined in terms of the matrix 
elements of the QCD currents between the vacuum and the one meson particle states 

(0 I j; I n') = -6-'e,m^F^ . (2.37) 



By using the KK decomposition and the currents (2.16) one finds, 

f 2 i'2 

^ (b' ± 6^) = ^ 

where -|- corresponds to the vector decay constant and — for the axial (obviously the only non 
zero matrix elements are between vector (axial) currents and vector (axial) mesons). 



T^K±b^) = §^b\, (2.38) 
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Finally the pion radius is defined in terms of the pion form factor as 



where 



with q^i = Pfi — p'u,- The pion form factor can be expressed as, 



K 



ft ' ^ m2-g2 



from which one finds that, 



K 



n=l " 



(2.39) 



(2.40) 



(2.41) 



(2.42) 



Actually a sum rule relates to Lg and /^r (see eq. (2.44)) so that (r^) = 12Lg/ f^. 



2.3.3 Sum Rules 

The discrete models satisfy several sum rules [6] (see [11] for the continuous case), all of which 



can be proven using the definitions given in the previous sections, the wave-equation (2.33) and 
the orthogonality of the eigenmodes. 

For example (generalized) Weinberg sum rules 



K p2 

n=l " 
K 



4Lc 



n=l 

f2 i-ir' ml^ F^^ = 0, 



(2.43) 



n=l 



and the following relations 



K 



n=l 
K 



n=l " 
K 



\ ^ FnjjB^ 

2 ^ ml 



(2.44) 



n=l 



hold independently of the number of sites. 
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Other sum rules however receive corrections for finite K. In particular vector meson dom- 
inance (VMD) is exact only in the continuous limit. From the fact that J^7r(0) = 1, form eq. 



(2.41 ) one has 



n=l 



To reproduce VMD fi should be significatively larger than As we will see this trend is 
reproduced by our fit of the data. 

The other sum rule which will be relevant in our phenomenological study is given by 



K r, A' + l 



This equation implies a deviation from the so called KSRF formula [23] , the phenomenological 
relation gp^^f^/m^ = 1/2, which is well verified experimentally. From the equation above it 
follows that this ratio is always less than 1/3. In comparison to continuous models, for finite 
number of sites there is an extra negative contribution to this relation which however quickly 
goes to zero as K oo. Both relations ( 2.45[2.46 ) hint to the fact that a realistic model of 



QCD will require at least a few sites to reproduce the data. 



3 Results 

The aim of this section is to quantify the agreement with the data of several continuous and 
discrete models. All the results are summarized in table [6] at the end of the paper. The 
main outcome is that many different models, including fiat space, have a good agreement with 
experimental results. Indeed, as we will see, only vague features of the metric are required to 
get fair predictions. 

In order to compare different models the first question to address is how to quantify their 
agreement with data. This is not obvious and involves a certain arbitrariness because the 
observables we compute are subject to a theoretical error that we can only estimate through 
power counting arguments. From this estimate we would expect typical relative errors in the 
10 to 30% range, larger than the experimental error for many observables. As in [7, 13] we 
proceed as follows: We restrict to the more precisely measured observables (namely, to those 
whose relative experimental error is less than 10%) for which the experimental error can be 
neglected in comparison with the theoretical one. We then define our error by the root mean 
square error estimate (RMSE). Denoting with Ol and Ol^^ the experimental and theoretical 
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values of each observable, we have 



RMSE 



1 ^ /C* - C' ' ^ 



E . (31) 



where No is the number of observables. The RMSE depends on the parameters of the model 
which we fix to their best-fit value by minimization; the minimum value of the RMSE is our 
measure of the error and is reported in the tables which follow for each different model. 

The statistical meaning of this procedure is the following: assuming that the theoretical 
error of each observable is AOl^ = ^Ocxp "^^ have 



RMSE = ^ 



1 ^f OL,-ok y , / iVo^ 



where d is the number of fit parameters. For a good quality fit Xmd ~ 1 so the RMSE {d will 
be significantly smaller than No) determines the size of the corrections necessary for agreement 
with the data. The assignment of the common relative error ^ is the main source of ambiguity in 
the procedure; there is obviously no rigorous way to establish whether a ^ error should be given 
to an observable or, say, to its square (9* . This translates in an ambiguity in the choice 



of the list of observables to be used in eq. (3.1). Our criterion was to use as observables the 



parameters of the corresponding matrix elements. We have checked that other choices, which 
of course produce slightly different results, leave the general conclusions unchanged especially 
for what concerns the comparison among different models. 

For each model we computed a total of 18 observables, 13 of which fulfill the "precision" 



criterion explained above and are used to compute the RMSE error of eq. (3.1 ). The 5 additional 
observables (shown in table |6| that we compute but do not use in computing the RMSE are: the 
Li, L2, L3 chiral Lagrangian parameters, which are not well enough measured, the p' mass and 
the Oi decay constant. We expect the p' (being it heavy) to be subject to larger theoretical errors 
than the other observables and this is why we exclude its mass from the RMSE. Concerning the 
Faj, we do not use the usually quoted experimental value of 144 MeV [24] because this value 
is extracted from the r decay in a model which is different from ours. Indeed, by computing 
the T — s> Stcu rate in our case (in which the branching ratio of ai is one and no other 

resonances contribute to the process) we find a slightly bigger value, around 165 MeV. Moreover, 
lattice calculations give a central value of Fa^ = 170 MeV [25]. All this suggests that Fa^ is not 
known precisely enough to be included in the RMSE. 
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3.1 5d models 



In the context of 5d geometries we present 4 representative models: flat space, AdS, Sakai- 
Sugimoto and for comparison a decreasing warp factor less steep than AdS, a{z) = ^JlJz. For 
a given geometry the model depends only on two free parameters: and L (as explained in 
section 2 we fix the gauge coupling of f/(l) and SU{2) to be equal). Even though for flat and 
AdS space analytic formulae can be used, in all the cases we compute the observables by first 



discretizing the model (through the dictionary (2.19) as discussed in section 2.2) and solving 
numerically. This method turns out to be very efficient numerically because the problem reduces 
to the diagonalization of a i^' x tri-diagonal matrix and allows to solve easily essentially any 
metric. Even for K = 500 (much beyond the required accuracy) the full fit can be performed 
in a few seconds. 

We now present the result for the most relevant cases: 



3.1.1 Flat Space 

The result of the global fit is reported in the table below. The main source of error of the model 
arises from the flat space relation rriai = 2mp. This effect is also enhanced by using square 
masses. 







mj 




fn 


Fp 




QpTTTT 






Qujp-K 




10^^10 


r^2 


Exp. 


.782 


.782 


1.22 


.087 


.15 


.14 


6.0 


.22 


.72 


15 


6.9 


-5.5 


12 


Th. 


.682 


.682 


1.42 


.081 


.12 


.12 


4.8 


.23 


.68 


13 


5.8 


-5.8 


11 


Dev. 


-23 


-23 


30 


-6 


-23 


-18 


-20 


3 


-6 


-11 


-15 


6 


-12 



Table 1: Flat space. The global error is 17%. All the dimensionfull quantities are in powers of 
GeV and the deviations are in percent. Fitted parameters: ~ 430 MeV and 5^2 ~ 7,5 MeV. 



3.1.2 AdS Space 

The global fit is. 

We see that the fit of AdS provides a better description of QCD data (particularly because 
of the right mass relations) even though fiat space could also be considered satisfactory at the 
level of precision we are working. However the superiority of AdS should not be overestimated. 
First of all our absolute measure of the agreement is subject to a certain arbitrariness. Our 
estimate of the error assumes that all the theorical errors are equal, and also depends of the 
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mj 




u 


Fp 




QpTTTT 








lO^Lg 


10=^Lio 


r^2 


Exp. 


.782 


.782 


1.22 


.087 


.15 


.14 


6.0 


.22 


.72 


15 


6.9 


-5.5 


12 


Th. 


.762 


.762 


1.22 


.082 


.16 


.16 


5.3 


.25 


.74 


15 


6.3 


-6.3 


11 


Dev. 


-4 


-4 


3 


-5 


6 


13 


-11 


13 


3 


3 


-9 


15 


-7 



Table 2: AdS space. The global error is 8.5%. All the dimensionfull quantities are in powers 
of GeV and the deviations are in percent. Fitted parameters: ~ 320 MeV and g^'^ ~ 5.3 MeV. 



choice of observables (for example square masses rather than mass). Moreover as we discuss 
in section [4], AdS requires a slighly lower cut-off than flat space as indicated by the value of 
in the fit. Having said this, the data do show a preference for a decreasing warp factor and 
AdS is already very close to the intrinsic error of the geometrical models (see below). In fact 
we showed by sampling a large number of geometries that the error is always larger than 7-8%. 
More in general about any gently globally decreasing metric will lead to an error less than 20%. 
For example a{z) = \/L/z gives a fit as good as AdS (see the summary table [6|. 

3.1.3 Sakai-Sugimoto Model 



For the case of the Sakai-Sugimoto model whose metric is given by (2.11) we find 





mp2 


mj 




u 


Fp 


F^ 




dpTT"/ 






lO^Lg 




r^2 


Exp. 


.782 


.782 


1.22 


.087 


.15 


.14 


6.0 


.22 


.72 


15 


6.9 


-5.5 


12 


Th. 


.782 


.782 


1.22 


.080 


.18 


.18 


5.5 


.26 


.79 


17 


6.6 


-6.6 


12 


Dev. 


-1 


-1 


-2 


-8 


22 


31 


-8 


20 


10 


12 


-4 


20 


2 



Table 3: The Sakai-Sugimoto model. The global error is 14%. All the dimensionfull quantities 
are in powers of GeV and the deviations are in percent. Fitted parameters: ^ 260 MeV and 
c/^2 ~ 5.4 MeV. 



It is interesting to see that the error is similar to fiat space. The main discrepancy is in this 
case the decay constant of the uj (here taking a different coupling constant for the U{1) would 
help to slightly fix this problem but would not improve the global error drastically). With 
respect to flat space this model reproduces more precisely the masses. 
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3.2 Few sites 



We now turn to purely 4 dimensional models. In comparison to the model proposed by Georgi 
in the context of the vector limit [3] we show that the addition of a few sites improves the 
agreement with QCD substantially. 

As we discuss in the next section, to reproduce QCD the five-dimensional theories considered 
have necessarily a low cut-off so that only a few resonances lie below the cut-off of the theory. 
Since our effective description describes reliably only the first 3 or so Kaluza-Klein modes, it 
makes sense to directly use a four dimensional model with few sites where only modes below 
the cut-off are included. The difference between this and an appropriate 5d model amounts to 
higher dimensional operators suppressed by the cut-off scale. For K sites {i.e. K resonances) 
the model has K-\-l parameters, the local cr- model decay constants and the gauge couplings (as 
in the continuous case the parameters of the C/(l) and SU{2) sector are taken to be identical). 

3.2.1 2-Sites 

Fitting the 3 parameters of the model we get the table below. The 19% error of this model is 
higher than the typical geometric model with decreasing metric. Note the poor predictions of 
Lq and of the pion radius. 







mj 




fn 


Fp 




QpTVTT 








10-^9 


lO'^Lio 


r^2 


Exp. 


.782 


.78^ 


1.22 


.087 


.15 


.14 


6.0 


.22 


.72 


15 


6.9 


-5.5 


12 


Th. 


.772 


.772 


1.12 


.085 


.12 


.12 


4.7 


.21 


.64 


11 


4.7 


-4.7 


8 


Dev. 


-3 


-3 


-23 


-3 


-17 


-12 


-21 


-2 


-10 


-25 


-32 


-15 


-35 



Table 4: Discrete model with K=2. The global error is 19%. All the dimensionfull quantities 
are in powers of GeV and the deviations are in percent. Fitted parameters: f\ k, 180 MeV, /2 ~ 120 
MeV and gi « 8.8. 



3.2.2 3-Sites 

In this case the number of free parameters is 4. The best fit corresponds to a 12% error. 

Therefore a model with 3-sites and nearest neighbor interactions provides an excellent fit 
of the QCD observables considered. For comparison, had we performed an analog fit in the 
Georgi's model [3] (clearly removing certain observables such as the ones associated to the Oi) 
the error of the fit would have been more than 30%. Note that the best fit drives the decay 
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mj 




U 


Fp 




QpTVTT 








lO^Lg 


10=^i^io 


r^2 


Exp. 


.782 


.782 


1.22 


.087 


.15 


.14 


6.0 


.22 


.72 


15 


6.9 


-5.5 


12 


Th. 


.772 


.772 


1.12 


.082 


.14 


.14 


5.0 


.23 


.70 


13 


5.3 


-5.3 


9 


Dev. 


-3 


-3 


-14 


-5 


-4 


3 


-17 


6 


-3 


-15 


-23 


-3 


-21 



Table 5: Discrete model with K=3. The global error is 12%. All the dimensionfull quantities 
are in powers of GeV and the deviations are in percent. Fitted parameters: /i 270 MeV, /2 ~ 130 
MeV, gi PS 7.3 and 52 ~ 10. 

constants in a region where they are hierarchical. This somewhat resembles the case of the AdS 
metric. In fact (see next section) the low value of the decay constant at the infrared worsens 
the calculability of the theory, but on the other hand a non generic choice of parameters is 
required to improve the fit, so there seems to be a balance between calculability and low error 
of the model. 

We also considered the case of 4 sites. In this case we find an error less than 10%, as low 
as the one for AdS (see the general table ([6])). We should note however that such agreement 
emerges at the price of rather large values for some of the g^s and again hierarchical values of 
the /'s. Increasing further the number of sites does not make the error decrease much further 
so we conclude that low energy QCD can be well described by a discrete model with 3-4 sites. 

3.3 The intrinsic error 

All the models that we have considered share common predictions, like for instance the relations 
between the Li coefficients of the O^p"^) chiral Lagrangian, the sum rules discussed in sect. 2.3.3 
and the equality of the p and u masses. These predictions are phenomenologically quite suc- 
cessfuU and partially explain why different models show a comparable and fair agreement with 
the data. Our model-independent relations, however, are not perfect and therefore produce 
an irreducible intrinsic error which is impossible to eliminate, even allowing for completely ar- 
bitrary geometries. We quantify this intrinsic error as 7-8%, which is the minimum error we 
obtain by scanning the parameter space of the discrete model with many {K = 100) sites. The 
intrinsic error is not far from the one of AdS or ^jL/z, meaning that these geometries have the 
non-trivial property of practically minimizing the error. 

It is natural to ask what fraction of the total error has to be really ascribed to the choice of 
the geometry and what part is due to the unavoidable intrinsic error. To answer this question 
we consider a restricted set of observables for which the tension due to the model-independent 
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predictions is eliminated. These are the masses of the p and of the ai, the decay constants 
of pion, p and uj and gu}pTv The fit gives now an error of 4.8% and 6.5% for AdS and 
respectively, while for the Sakai-Sugimoto model the error practically does not change and for 
flat space it even increases by about 1% due to the bad mass relations. The scan of the discrete 
model parameter space gives now a minimum error which is essentially zero, showing that indeed 
we have completely eliminated the intrinsic error with this new choice of the observables. These 
results have been obtained by using more precise experimental values and allowing for different 
couplings in the U{1) and SU{2) sectors, i.e. by fitting the parameter a of eq. (2.8) instead 
of fixing it to one. At this level of precision small deviations of a from one have indeed a 
significant impact. 

Some of the relations which produce the intrinsic error are obvious: the mass of the uj is the 
same as the mass of the p, Lg and Lio are fixed by r^, and and the couplings gp.„^ and g^^^,^ 
are expressed in terms of (^p^^r and /^r (see sect. 2). Another relation, already discussed below 
eq. ( |2.46[ ), is that 



rUp^ 



< 1 , (3.3) 



independently on the warp factor, while the experimental value is about 1.3. Also, the sum 
in eq. (2.42) which gives the radius of the pion is essentially always saturated by the first 
resonance. We therefore find the relation 

%^J^ = l, (3.4) 

but the experimental value is so close (1.02) that even the small deviations due to the rest of 
the sum in eq. (2.42) allow the experimental result to be exactly reproduced. 

It is important to remark that the model-independent predictions discussed in this section, 
which lead to the minimal intrinsic error, are only valid at the tree-level and will be modified 
by higher order corrections. A non-zero value of Lg + Lio, for instance, will be generated by 
loop corrections and its size is expected, by simple power counting arguments, to be consistent 
with the experimental one. The vanishing of Lg + Lio at tree-level is related with the tree-level 
absence of the oi 717 decay, which is another feature of our model. This relation is lost, 
however, at the loop level and indeed even in the HLS case, in which the Oi meson is not 
present, a non-vanishing Lg + Lio is generated [2]. It would be interesting to etablish whether 
radiative corrections induce a realistic ai —>■ TC'-f decay in the models we have considered. 
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4 Calculability 



One important feature of extra-dimensional descriptions of QCD is provided by the calcula- 
bility of these models. We now wish to consider this aspect in detail. The arguments below 
become good estimates in the case of a large number of colors. We consider the case of a 
flat geometry where the power counting is transparent. Naively, neglecting boundary effects, 
the five dimensional theory has a regime of validity up to the scale A ^ IQn'^/g'^. Given the 
relations, 

9 



9iL 



= ^ (4.1) 

and the fact that 47r/^ ^ y/Ncirip it follows that l/t^g is linear in A^^- In absence of other 
interactions this would imply that the number of weakly coupled resonances below the cut-off 
scale is also roughly Nc. 

Due to the presence of the Chern-Simons term, necessary to reproduce the anomalies of the 
theory, the number of weakly coupled resonances is however smaller for A''^ large, see also [13]. 
This is due to the fact that the coefficient of the Chern-Simons term scales with Nc. This 
interaction is more relevant than the pure non-abelian gauge interaction at large Nc and lowers 
the maximum possible cut-off of the theory. Considering the growth with energy of scattering 
amplitudes of gauge bosons in 5D the new cut-off can be estimated to be 



K 



2 



3 167r2 



where k is a numerical factor of the order few for the most relevant interaction (see eq. (2.8). 
This means that the maximum number of weakly coupled resonances that the higher dimen- 
sional theory can possibly describe scales as 

A^maxOcArV3. (4.3) 

What is important for the calculability of the 5d theory is that N^ax still grows, allowing to 
have a weakly coupled deformation of the theory by taking Nc large. Despite the fact that the 
Chern-Simons lowers the cut-off, the couplings of the mesons respect the correct scaling of large 
Nc QCD. This follows from the fact that since l/gf is proportional to A'c we can factor out the 



^Independently of anomalies in the gravitational dual of theories such as = 4 Super- Yang- Mills one finds 

2 

M^/rUs oc Nc where is the string scale. It is therefore reasonable to expect in general that this cutoff 
corresponds to a physical scale. 
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number of colors in the five-dimensional action. As a consequence interactions of the KK modes 
are suppressed at large Nc and the appropriate countings are reproduced. Using the previous 
rules for the case of A^^c = 3 considered in this paper we note that the Chern-Simons does not 
lower cut-off in this case and 3 resonances are weakly coupled in our effective description. 

It is interesting to see how the above results can be understood in the deconstructed version 
of the theory. The basic equations in this case are given by 

f2 _ J_ 

- X' 

= ^, (4.4) 

where K is the number of sites. It follows that IGn'^K ^ g'^N^. Since g < 4tt for the gauge 
theory to make sense, this already tells us that the 4d description cannot contain more than 
0{Nc) states that are weakly coupled. However, as in the 5d case, this upper bound cannot 
be saturated due to the interactions associated with the anomalies. The cut-off can be com- 
puted locally by looking at the scale where the individual cr— models become strongly coupled. 
In the discretized version the Chern-Simons term decomposes into gauged WZW terms with 
coefficients Nc at each site (see appendix). The most relevant interaction arises from the term 

which becomes strongly coupled at the scale A ^ GAn^k f /{Ncg"^). Requiring that the heaviest 
state lies below the cut-off one finds that the number of weakly coupled resonances is maximized 
by taking 

An 

g^—. (4.6) 



This implies, as in the continuous version, the relation (4.3). Note that contrary to the pure 
gauge case reproducing the 5d cut-off demands that the local gauge coupling goes to zero. 

Therefore the geometrical picture and the deconstructed version with nearest neighbor inter- 

1 

actions and a number of sites ~ Nc are indeed equivalent as effective field theories. 

Similar considerations can be repeated in the case of a generic geometry /deconstructed 
model. An interesting fact is that, working with rup and fixed, fiat space actually maximizes 
the cut-off defined by the 5d gauge coupling. For example to reproduce the experimental values 
of and nip AdS requires g"^ which is 30% larger than fiat space. In fact the result of our 
fit is strongly correlated to the dependence of the formal cut-off of the theory. Of course such 
statements should be taken with great care as there is no parametric separation between the 
cut-off in fiat and AdS space. However for the case of interest, Nc = 3, where the calculabilty 
of the model is inevitably limited, such a numerical difference has an impact. 
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5 Discussion 



In this paper wc compared several different holographic models of QCD to quantify their agree- 
ment with data. It is rather remarkable that, compatibly with the corrections expected in these 
type of models, essentially any smooth decreasing metric and even fiat space provide a good 
fit of the low energy mesonic observables. This suggests that there is no special virtue in the 
choice of the AdS metric originally inspired by the AdS/CFT correspondence. Nevertheless it 
is interesting to note that, in the class of models we considered, AdS almost minimizes the error 
while flat space and Sakai-Sugimoto models give a similar error which is higher than AdS. 

Perhaps the most interesting outcome of our analysis is that the successful results of the 
holographic models of QCD arc already reproduced by a four dimensional model with 0{3) 
sites. This is to be expected because all the flve-dimensional models are characterized by a low 
cut-off flxed by the experimental values of rUp and fj^. As a consequence only 0{3) resonances 
can be reliably described within the effective theory and a macroscopic notion of geometry is 
lost. In this sense it might be useful to consider holographic theories of QCD as discrete theories 
containing only the modes lying below the cut-off of the theory. These are generahzations of the 
model proposed by Georgi in [3] . What is crucial phenomenologically is the structure of nearest 
neighbor interactions of the discrete theory/extra-dimension which also ensures an enhanced 
calculabihty of the theory. 

Among possible extensions of our analysis the most important would be the inclusion of 
the baryon sector of the theory. As shown in [13] the baryons arise in the 5d picture as 
topological solitons analogous to 4d skyrmions. The extra-dimension allows again to make 
such configurations calculable contrary to their 4d cousins. The same will be true for any 
metric and discrete model considered here, even though explicit solutions might not be easy to 
find. As in [13], it would be interesting to consider the impact of the baryon observables on the 
global fit. 

We do not expect one model to be strongly preferred. Due to the small value of Nc in 
QCD the calculability of the theory is always limited so many different models will always 
give comparable predictions. In this regard the AdS/QCD like models considered in this paper 
make sharp predictions for a world with a large number of colors already at tree-level, so if such 
data were available from lattice computations, it would allow to test whether such approach 
has some truth or should be rejected, and which geometry is favored. Alternatively one could 
try to compute next-to-lcading order corrections, paying the price of adding new independent 
parameters, to check whether the intrinsic error of these models is reduced and if a preference 
for some geometry emerges. 
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A Anomalies and the Wess-Zumino-Witten term 



In this appendix we show how the QCD anomalies can be incorporated in the discrete model 
by the introduction of a gauged WZW at each link. Due to anomalies the effective action of 
large- Ac QCD must transform as 



6S 



247r2 



(A.l) 



where the 4-form uj\{cx.,A), whose exact definition is irrelevant for the moment, represents the 



QCD anomaly and will be given in eq. (A. 11). Infinitesimal local chiral symmetry transfor- 
mations, which correspond to L = Li and R = Hr+i in eq. (2.12), are parametrized by cx.lji 
as L ~ 1 + iai, and similary for R. Notice that in writing eq. (A.l) we have ignored the 
U{1)aSU{Nc)c-SU{Nc)c anomaly, which is the one responsible for the rj' mass. This is justi- 
fied in the large- Ac limit where the latter anomaly is subleading [26] and the t]' is massless as 
it will be in our model. Needless to say, the anomaly is only defined up to the variation of local 



counterterms, so that we had to make a choice to write eq. (A.l). With our choice the anomaly 



is entirely in the U{1) part {i.e., it vanishes for SU{2) transformations, see eq. (A. 11)) and 



assumes the "Left minus Right" form. Having made a choice does not imply, however, a loss 
of generality: any other form of the anomaly could be obtained by adding local counterterms 
to the Lagrangian we will derive assuming eq. ( A.l[ ). 

According to the previous discussion, an "anomalous" non gauge-invariant term, with vari- 



ation given by eq. (A.l), has to be added to the action; we will denote the latter as Sa- In 



order to construct Sa let us assume that each of the K + 1 a-models is anomalous, exactly as 
it would be if it represented the effective description of a large-Ac QCD-like theory. If this is 
the case the action, prior to the gauging, contains a gauged WZW term for each link 



K+1 



WZW 



(A.2) 



A;=l 



^If we had, for instance, to put the anomaly in the standard form in which the vector subgroup is anomaly- 



free we would just add the Bardeen counterterm of eq. (A. 15 1 
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where Twzw will be defined in eq. (A. 16). The WZW terms, according to eq. (A. 6), produce 
an anomaly 



K+l K+1 
k=l k=l 



(A.3) 



Notice that the variation (A.3) vanishes for a global transformation and eq. (A. 2) is still invari- 



ant under the entire global group of eq. (2.12). 



We must now, as we did before with the Lagrangian (2.13), gauge the vector combination 
of the 1'^'''^ and sources. To this end we replace l'^ with A'^"^ and with A'^ in the action 



(A. 2) and obtain 



K+l 



Sa = ^ r^/zvy [5]^, a'' \A''] 



(A.4) 



k=l 



Under the local groups that we would like to gauge and under the ungauged "end-points" 
groups which correspond to the chiral QCD ones the action 5*^ transforms as 



6S, 
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[j^Pl{a'-\A'-')-uJl{cx\A')] = ^ / K(«L,l)-^l(«^,r)] .(A.5) 

k=i 



The second equality follows from the identifications cx^ 



K+l 



1 and 



A^^^ = r and from the fact that the anomaly of each group, except the first and the last ones, 
receives equal and opposite contributions from the links to its left and to its right. This ensures 
that the anomaly is cancelled for the vector groups, which can therefore be safely gauged, and 



that Sa has the correct anomaly required by eq. (A.l). Moreover, thanks to our construction. 



Sa is in the nearest-neighbor form of eq. (2.18) like the rest of the action. Notice that for the 



cancellation of the vector anomaly it has been crucial that we chose the anomalies of all the 
cT-models to be equal {i.e., to correspond to the same number of colors Nc) and each of them 
to be in the "Left minus Right" form. It is easy to realize that if we had started from any other 
form of the anomaly, as long as it was of the "Left minus Right" type, we would have obtained 
the same Sa action. 



We now need, in order to write the explicit form of Sa of eq. (2.17), to work out the Twzw 
action. This is the subject of the next subsection. 



Gauged WZW for U(2) 

The general expression for the gauged WZW term in the case of Nj flavors, i.e. of U{Nf) 
or SU{Nf) chiral groups, is well-known and can be found in [27]. In the two flavors case 



26 



that we are considering, however, the general formula can be strongly simplified and, up to 
gauge-invariant operators, the WZW assumes a much simpler form. It is more convenient and 
more interesting, instead of starting from the general formula and simplifying it, to derive from 
scratch the WZW term for two-flavor QCD and this is what we will do in the following. We 
will use a form notation in which A = —iA^dx'^, where A denotes here any Left or Right gauge 
field, and the field strength is a 2-form F = dA + A^. 

Our starting point is the chiral QCD anomaly which, by a standard text-book calculation 
can be written as 



(A.6) 



where ST denotes the variation under local chiral transformations cxl,r of the effective action 
for the U{2)l x f/(2)/j sources 1 and r. The uj\ 4-form is defined as 



OL 



d (^AdA + ^ A^^ ^ , 



(A.7) 



and eq. (A.6) provides the standard "symmetric" form of the anomaly. It is customary to define 
the Chern-Simons 5-form uo^ by regarding the ordinary 4d space-time as the boundary of a 
fictitious 5d and extending the sources to 5d gauge fields in this space. It is 



^5(A) 



A (dA)' + -AMA + -A^ 
^ ' 2 5 



(A.8) 



and the variation of uj^ is the exterior derivative of the anomaly: ^^oj^ = dujl- 

All the equations above hold in general, for any number of fiavors, but an important simpli- 
fication occurs if we specialize to the case of Nf = 2 in which A = 1/2A + A is an U{2) gauge 
field. The CS becomes 

1 



^5(A) = z^i(dA)2 + z^l(F2) + zd 



A( AF 



4 



-A' 



= lJ.(A) + dX. 



(A.9) 



This means that, by adding the local counterterm X to the effective action {i.e, by a change 

.nomaly 

6T [1, r] 



of the regulator), the anomaly (A.6)) can be put in the form 
where the 4-form 



(A. 10) 



^l(a,A) = \a(dAy + la{F') , (A.ll) 
is defined, up to an irrelevant total differential which cancels under integration, from the relation 

5aUJ5 = dujl. (A. 12) 
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In the alternative form (A.ll ), the anomaly is entirely in the U{1) part and the effective action 



r is SU{2)l X SU (2) n-mYanant. 

In the low-energy effective action for the Goldstone bosons U{2) matrix S the anomaly 



(A. 10) is reproduced by a suitably designed gauged WZW term, T^zw P,l,r], which is such 



that its variation respects eq. (A.ll). Following ref. [28], a term with this property is 



WZW 



[S,l,r] 
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I {[^5(1) -^5(r)] - [^5(l('^'^)-^5(r) + dfi4(l('^'\r)]} , (A.13) 



where 1*^^^^ denotes the Left source gauge-rotated by the inverse Goldstone matrix S"!". Given 
that S transforms as S — g^'SgR^ , l*-^^-* transforms in the same way as the r source does. 



The way eq. (A.13) works is not difficult to explain: the first square bracket alone already 



gives the anomaly but is not a closed form so that it depends on the unphysical extention of 
the sources (and of the Goldstones) in the extra dimension. The second square bracket must 
be there to give a closed form when summed with the first one, and must be gauge-invariant. 
It clearly accomplishes the first task as duJ^ = du^ = (F^) but for the second one we have 
to choose a suitable "Bardeen counterterm" 4-form B4. This must be such that under vector 
transformations (remember that l'-^^-' transforms as r) uJ^i}) — cJ5(r) + di?4(l,r) is invariant. 
This means that 



5y54(l,r) 



uj](a., r 



(A. 14) 



Given that the anomaly is enterely in the U{1) part, B4 must be invariant under SU{2) vector 
transformations and moreover, due to parity invariance, it must be odd under the 1 ^ r 
operation. Also, it must contain only one derivative as its gauge variation must contain two. 
The only two independent terms with those feature which contain SU{2) fields and are not 
gauge invariant are iv{a{L + R)) and iv (a^) where v = (1 + r)/2 and a = (1 — r)/2. The 
coefficients are extracted by direct calculation and the result is 

54(1, r) = - J (Tr{dT+ df)) + {T+ f) ((/ - r) [3 (L + R) - {I - rf] ) . (A. 15) 



Applying now eq. (|A.13|) and after some simple algebra we can write 

Tr{di + dr) - m ({dTy 



WZW 



[S,l,r] 



967r2 



(df) 



+ 



(?+ f) /3 (^LSDS"^ - + 



d/df ) -3m(L2 
3^ 



R') 
(A.16) 



where we have separated S = SS in the U{1) and SU{2) component and parametrized S as 
S = exp zll. Notice that, as known, no ungauged WZW term is present in the U{2) case and, 



consequently, T\yzw in eq. (A.16) is enterely written as an ordinary 4d integral. 
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The WZW term, in the usual NDA counting of the chiral Lagrangian in which gauge fields 
and derivatives both count as 0{p), is an O(p^) term and its coefficient is perfectly consistent 
with the NDA rule. |^ At the same order, several gauge-invariant operators appear, which are 
listed in general in ref. [29] for the U{3) chiral group. The WZW term, therefore, is in any 
sense unique given that we could have added to eq. (A. 16) any combination of these 0{p'^) 
gauge-invariant operators. 



*One has to remember that physical cut of the cutoff of the cr-model, i.e. the scale at which the p meson 
enters, is A^sb = 47r/^/\/]Vc at large- A^c, and not 47r/^. 
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Table 6: All the dimensionfull quantities are in powers of GeV and the deviations are in percent. 

The fitted parameters (MeV): for flat space, ~ 430 and g^^ ~ 7.6; for AdS space, ~ 320 and g^^ ~ 5.3; for the Sakai- 
Sugimoto model, L^^ « 260 and g^'^ ^ 5.4; for -y^L/z, 370 and g^'^ 6.3; 

for K = 2, fi^ 180, /2 ^ 120 and gi ^ 8.8; for if = 3, /i « 270, /2 « 130, gi ^ 7.3 and g2 ~ 10; for if = 4, /i 460, /2 fa 150, 
fs 140, 51 5.3 and 52 ~ H- 
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